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Abstract 

Quantum vorticity occurs in superfluidity, which arises from a spatial variation 
of the quantum phase. As such, it can occur in diverse systems over a wide range 
of scales, from the electroweak sector and QCD of the standard model of particle 
theory, through the everyday world, to the cosmos. I review the observable man¬ 
ifestations, and their unified description in terms of an order parameter that is a 
complex scalar field. 


1 Overview 

Quantum vorticity is a manifestation of superfluidity, which arises from a spatial varia¬ 
tion of the quantum phase. This can occur in different types of physical systems, but a 
unified description can be given in terms of an order parameter that is a complex scalar 
field; 

0(x) = (1) 

Its existence signals the spontaneous breaking of a t/(l) gauge symmetry. 

In atomic systems with Bose-Einstein condensation HI, the broken gauge symme¬ 
try is global. In this case, represents the condensate wave function, which can be 
described by the Gross-Pitaevskii equation, a nonlinear Schrodinger equation (NLSE): 

+ ( 2 ) 

where m is the mass of the atoms, g an interaction parameter, and /uq is the chemical 
potential. 

In superconductivity ||2l, the broken gauge symmetry is local, and <p obeys a non¬ 
linear Schrodinger equation with gauge coupling: 

tP' / ifiq , d(p 

- ^ ( V - 1 (p + ig^ (p- juo) (p = in— (3) 
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where q — 2e is the charge of Cooper pairs, and A is a magnetic field. This is called 
the Ginsburg-Landau (GL) equation. 

In relativistic systems Q, these equations generalize to a nonlinear Klein-Gordon 
(NLKG) equation : 

□0 + V'cf>^J (4) 

where □ is the d’Alembertian in curved spacetime, which may contain gauge couplings, 
and J denotes a possible external current. The self-interaction potential may be chosen 
to be a phenomenological potential; 

y = ^ ( 0 V - Faf + Vo ( 5 ) 

Thus 

V'^A(f<p-Fo) (6) 

where a prime denotes differentiation with respect to 0*0. Physically, 0 may corre¬ 
spond to a component of the Higgs field in the standard model, the order parameter 
associated with chiral symmetry breaking in QCD, or Higgs-like fields in grand-unified 
or supersymmetric theories. 

All these systems exhibit superfluidity, with superfluid velocity given by Q 


^ kVct (7) 

where k - -c (dcrldty^ . In the relativistic case, the requirement |Vj| < c is guaranteed, 
when is time-like, i.e. ,d^crd^cr < 0. In the nonrelativistic limit, cr - {mc^ lfi^t + 
v,and K tijm. 

The equations of motions can be rewritten in terms of F and cr, and the equation 
for cr yields a equation for Vj similar to the classical Euler equation of hydrodynamics, 
with quantum corrections. However, for many purposes it is simpler to stay with the 
original forms, particularly numerical analysis. 

Quantum vorticity arises from a quantization of the circulation, by virtue of the 
continuity of 0: 


dx ■ V cr - Inn (n = 0,+1,+2,...) (8) 

c 

where C is a closed loop. If n 4 0, the loop encircles a vortex line on which 0 = 0. As 
illustrated in in Fig.l, the field rises from zero at the line to an asymptotic value, with 
a characteristic healing length The vortex line can thus be thought of as a tube of 
radius ^ in which symmetry is restored, with expulsion of of the order parameter. The 
superfluid velocity Vj drops off with distance r from the line center like r '. Since the 
energy density of the superfluid flow is proportional to Vj, the energy per unit length of 
of a single straight vortex line diverges with the radius of the container. The existence 
of vorticity depends on a vacuum field, hence on the nonlinearity; it is absent in the 
usual Schrodiger equation, or Klein-Gordon equation, where the field goes to zero at 
infinity. 

A vortex line cannot terminate, except on boundaries, or upon itself, forming a 
closed curve. The simplest example of the latter is vortex ring of radius R, as illustrated 
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Figure 1: The vortex tube. The superfluid velocity decreases inversely with distance 
from the central vortex line. 


in Fig.2. The ring moves with a translational velocity normal to the plane of the ring v ~ 
' Inf? , with energy E ~ R\nR. Thus, roughly speaking, E ~ v^'.A vortex line of 
arbitrary shape will move in such a manner that, at any point of the line, the local 
velocity is normal to the tangent circle, with velocity inversely proportional to the local 
radius of curvature. Thus, the line will generally execute a self-driven writhing motion, 
as illustrated in Fig.2. 

When two vortex lines cross, they will reconnect, as illustrated in Fig.3. In the 
final configuration, there appear two cusps, which spring away from each other at great 
speed, because of the smallness of the radii of curvature. A reconnection thus creates 
two jets, which represents an efficient way to convert potential energy into kinetic en¬ 
ergy in a short time. The reconnection of magnetic flux lines in the sun’s corona is 
believed to create jets of solar flares, as shown in Fig.3. 

Vortex rings can be created in a superfluid with a heat source. They would ex¬ 
pand, intersect and reconnect, and reach a steady-state vortex tangle called quantum 
turbulence. Fig.4 Q shows a computer simulation of the process. The vortex tangle 
is a new geometrical object, with a fractal of dimension 1.6 Q. When the heat source 
is removed, the vortex tangle will decay. The dynamics of growth and decay can be 
described phenomenologically by Vinen’s equation Q 



(9) 


where ( is the average vortex-line density i (length per unit volume), and A and B are 
phenomenological parameters, with A governing the growth, and B the decay of the 
vortex tangle. 

In the following, we briefly survey quantum vorticity in diverse physical systems. 
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Vortex ring 



Figure 2: Left panel: The vortex ring has a translational velocity v approximately 
inversely proportional to its radius R. Right panel: A general vortex line has local 
velocity normal to the tangential vortex ring, with magnitude inversely proportional to 
the radius of curvature. 



Figure 3: Left panel: Reconnection of vortex lines creates two cusps that spring away 
from each other at high velocity, creating high-velocity jets of flow. Right panel: Jets 
in magnetic reconnections in the sun’s corona creates solar flares. 
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Figure 4; Emergence of quantum turbulence in computer simulation. Numbers under 
each picture indicate the number of reconnections. From Ref. 0. 


2 Liquid helium and cold trapped atomic gases 


Quantized vortices in superfluid helium and atomic gases can be described by the NLSE 
•HI- with superfluid density given by n = and healing length given by 




/to 


( 10 ) 


By shooting alpha particles or electrons into superfluid helium, vortex rings can 
be created, which trap the projectiles |8l. The velocity-energy curve of the composite 
object can be obtained by dragging it cross the liquid against an electric held. The 
results fit that of a quantized vortex ring, as shown in Fig.5. 

In superfluid helium and cold condensed atomic gases, vortices can be created by 
rotating the container. When the angular velocity of the container exceeds a critical 
value, the superfluid responds by developing vortex lines parallel to the axis of rotation 
El nni. Experimental results in a trapped atomic gas are shown in Fig.6, showing the 
development of a vortex lattice as the angular frequency increases IfTTI . 

Quantum turbulence iflBl . as well as vortex reconnections ifT^ . have been studied 
in superfluid helium. The velocity distribution in quantum turbulence is found to have 
a power-law tail as shown in Fig.7. This is quite different from the Gaussian 
distribution of classical turbulence, and is due to the occurrence of large velocities 
in vortex reconnections. Computer simulations using the NLSE have reproduced this 
distribution KH. 

Vortex lines in liquid helium have been made visible by coating them with nano¬ 
sized metallic dust, as shown in Fig.8 na. The metallic dust adheres to the surface of 
the vortex tube because of the Bernoulli pressure arising from a decrease of superfluid 
velocity away from the the core. 

The computer simulations in Fig. 9 m demonstrates that a Bose condensate far 
from equilibrium creates quantum turbulence, and relaxes through its decay. 
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Figure 5; Velocity-energy curve of vortex ring created by ions shot into liquid helium. 
The ion (alpha particle or electron) gets trapped in the ring it created. The theory curve 
is that for a vortex with one unit of quantized circulation. 



Figure 6: Lattice of vortices in rotating cold trapped atomic gas, at increasing angular 
velocity of rotation. From Ref. m. 



Figure 7: Velocity distribution in quantum turbulence in superfluid helium (data 
points), as compared with the Gaussian distribution of classical turbulence (solid 
curves). The former has a tail due to the jets of large velocities accompanying 
vortex reconnections, which are essential for the maintenance of quantum turbulence. 
From Ref. ifTSll . 
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Figure 8: Vortices in superfluid helium made visible by coating the surface of the vortex 
with metallic dust of nano size. The adhesion is due to Bernoulli pressure. From Ref. 

ca. 


t = 1000 t — 2000 t = 3000 






Figure 9: Computer simulation of the evolution of a Bose condensate far from equi¬ 
librium, showing that it equilibrates through the creation and decay of quantum turbu¬ 
lence. These pictures show the decay sequence. From Ref. lfT^ . 
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3 Superconductivity 


In superconductivity, as described by the GL equation there is another length be¬ 
sides the healing length the penetration depth d. This arises from the conserved 
current ^ 

( 11 ) 

Zmi me 


When this is substituted into Maxwell’s equation V x V x 
Coulomb gauge V ■ A = 0, 


Annq^ 


A = 0 


A 


^J, we obtain, in 


( 12 ) 


for uniform 0*0 = n. The second term in the brackets is the photon’s squared mass 
inside a superconductor. An external magnetic field can penetrate a superconductor 
only to a depth d corresponding to the inverse mass: 





(13) 


This is the Meissner effect, a simple example of the Higgs mechanism, i.e., spontaneous 
breaking of local gauge invariance gives mass to the gauge particle. 

The presence of two characteristic lengths leads two types of superconductors: If 
^ > d, we have type I behavior: the magnetic field is completely repelled from the 
interior. If ^ < d we have type II behavior: the magnetic field penetrates the body 
in vortex tubes containing quantized magnetic flux heje, arranged in a regular lattice 
know as the Abrikosov lattice. We can understand these behavior by referring to Fig. 10, 
which depicts the interface between magnetic field and order parameter. In the type I 
case, there is little overlap between magnetic field and order parameter. If a flux tube 
were formed, the system can lower the energy by expelling it outside, and filling in 
the hole in the order parameter. For Type II, an overlap between magnetic field and 
order parameter lowers the energy, creating a negative surface tension. Thus, flux 
tubes are formed, penetrating the superconducting body, and maintained by solenoidal 
supercurrents. If the medium is infinite, the flux cannot be expelled, and finite-energy 
configurations are flux tube of finite length, terminated by magnetic monopoles. In this 
regard, see Fig. 11 for a comparison between the magnetic field due to two monopoles 
in vacuum and in a superconductor. If the tube is cut, one does not get free monopoles, 
but more flux tubes terminating in N and S. In this sense, a superconductor is a medium 
of magnetic-monopole confinement. 


4 “Vorticons” in the Higgs field 

In the electroweak sector of the standard model, the local gauge symmetry is sponta¬ 
neously broken by of the multi-component Higgs field, which gives mass to the gauge 
bosons W± and Zq, while leaving the photon massless. We can envisage a ’’vorticon” 
El, a microscopic vortex ring in the Higgs field, in the shape of a donut, containing a 
gauge boson that is massless, as illustrated in Fig. 12. It is a microscopic donut-shaped 
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Figure 10; Upper panels show the behavior of the magnetic and order parameter at an 
interface, for type I and type II superconductors, lower panels illustrate the expulsion 
of magnet flux in type I superconductors, and the formation of a quantized-flux lattice 
in type II case. 



Figure 11: The magnetic field of a magnetic dipole in empty space (left) and in a 
superconductor (right). In the latter the magnetic flux is squeezed into a flux tube by 
circulation supercurrents. The situation suggests that a superconductor would exhibit 
magnetic monopole confinement. 
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Figure 12: A Zq vector boson created in a high-energy collision induces superfluid flow 
in the background Higgs field, creating, and is trapped by, of a vortex ring, resulting in 
a unstable particle dubbed the “vorticon” , with a mass of approximately 3 TeV, and a 
lifetime of the order of 3xl0“^^s. 


waveguide, which could be created when an energetic Zo,say, tears through the Higgs 
field, much as a projectile can create and be trapped by a vortex ring in liquid helium. 

The mass of a vorticon can be estimated by constructing normal modes of the Zo 
field inside a torus cut out from the background Higgs field, and minimizing the energy 
of the lowest mode El. The standard-model Hamiltonian in the Zo sector is, with 


h - c - I, 



i(B2+E2) + |(V-/^Z)0|2 + F(0) 


(14) 


Here, cp is the Higgs field, Z the vector potential in Coulomb gauge (V ■ Z = 0), and 
B = V X Z, E = -(9Z/clf.The Higgs potential is given by 


with 


(15) 


Fq = 174 GeV 

y/lAFo — niH = 125 GeV 


(16) 


where mu is the Higgs mass. Thus A = 0.256. The gauge coupling constant q and Zo 
mass are given by 


q 

mz 


e 

sin 26]fi/ 
gFp 

V2 cos 6w 


(17) 

(18) 


where e is given through the fine-structure constant e^/A-n x 1/137, and 6w is the 
Weinberg angle given through sin^ 6w ~ 1/4. 

There are two types of vorticons: magnetic and electric, with the magnetic (electric) 
field pointing along the toroidal direction. As in electromagnetic wave guides, there is 
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no completely transverse mode. The masses are found to be 


^mag 

mz 

^elec 

niz 


35.3 + 6.42 



1.03 



27.7 + 5.66 


\2/3 

mH \ 

mz) 


+ 0.504 



( 19 ) 


With the experimental values mz -91 GeV and mn = 125 GeV, we have 

M,„ag = 3.47 TeV 

Melee = 2.88 TeV (20) 

The size of these vorticons are of order ~ 10 cm. They are unstable, with 
lifetimes of the order the Zq lifetime 3 x lO^^^s. 


5 QCD strings 

A phenomenological description of quark confinement in QCD can be modeled after 
the magnetic monopole confinement in a superconductor. The difference is that quarks 
generate Yang-Mills color-electric fields instead of magnetic fields, so the QCD vac¬ 
uum confines electric flux instead of magnetic flux. (In the nonlinear Yang-Mills gauge 
theory, there is no electric-magnetic duality as in Maxwell theory.) A meson consists 
of a quark and antiquark pair, belonging respectively to the representations 3 and 3 
of color SU(3), connected by a flux tube. This can be represented by the picture in 
Fig. 11., if we substitute 3 for N, 3 for 5. The order parameter corresponds to a con¬ 
densate of quark-antiquark pairs, and the electric flux tube is maintained by solenoidal 
color-magnetic currents. A baryon made of three quarks can be represented by a flux 
tube with one quark at one end, at two at the other. The two-quark system has color 
representations 3x3 =3-1-6, and we can use the 3 irreducible combination. Alter¬ 
natively, the baryon could be represented by Y-shaped flux tubes each terminated by 
3. 

One can approximate the flux tube by a massless string terminated at both ends by 
quarks, with the string rotating about a perpendicular axis, under tension Tq lUSl, as 
illustrated in the upper panel of Fig. 13. In the limit of small quark masses, one obtains 
a relation between between J and M: 

J = ao + a'M^ (21) 


where ao is a model parameter, and 


/ 

a 


1 

ItiTq 


( 22 ) 


As shown in Fig. 13, J-M^ plots of meson and baryons do yield straight lines known as 
Regge trajectories, with a universal slope a' ^ 0.9 GeV'^. This corresponds to a string 
tension of approximately 16 tons. 
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Figure 13; Hadrons modelled as strings terminated in quarks. The strings are idealiza¬ 
tion of vortex tubes containing color-electric flux. Their rotational angular momentum 
J and squared mass bear a linear relation that agrees with the observe Regge tra¬ 
jectory of particles, and show the string has a tension of about 16 tons. 


It would be desirable to set up an analog of the GL equation for QCD, but the sit¬ 
uation here is more complicated than superconductivity, because of the intrinsic non¬ 
linearity of Yang-Mills gauge fields. Besides quark confinement, there is chiral sym¬ 
metry breaking. The condensate could involve topological object such as monopoles. 
We refer the reader to CD, lEol, ED, for works related to that goal. 

6 Cosmology 

The Higgs field of particle theory has found experiment support with the discovery of 
the Higgs boson Il22ll . Il23l . Grand-unified or supersymmetric theories call for more 
Higgs-like fields. Any long-range order, such as the phase coherence of these fields, 
will persist to infinite distances in spatial dimension greater than two. In two or fewer 
dimensions, the long-range order will be destroyed by the long-wavelength fluctuations 
of the Goldstone modes ED- The Higgs field, and other possible Higgs-like fields, thus 
make the whole universe a superposition of different types of superfluid, (like a mixture 
of "^He and ^He below lO^^K.) We discuss the manifestations of superfluidity in terms 
of a generic complex scalar field, in three different epochs of the universe: the big-bang 
era, the CMB (cosmic microwave background) era, and the present. 

The scalar field presumably emerges during the big bang, but how it does this de¬ 
pends on model. A quantum field needs a high-momentum cutoff A, which should be 
infinite at the big bang. If the potential emerges from zero at that moment, it must be 
asymptotically free. This rules out all polynomial forms, and admits only the Halpem- 
Huang potential E5l . with exponential behavior for large fields. A big-bang model 
has been constructed based on this potential ED E3, with uniform Robertson-Walker 
metric whose scale is tied to the cutoff: a - A“'. Solving Einstein equation numeri- 
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Remnant vorticity 


Matter created thi'ough 
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Figure 14; The inflation era; creation of matter all matter the vortex reconnections in 
quantum turbulence. Ref. El. 


cally, as an initial-value problem, shows that the length scale expands like 



(23) 


where cq and p < I depend on model parameters and initial data. The result is equiva¬ 
lent to having a cosmological constant that decays in time according to a power law. If 
this behavior persists to later times, it could explain dark energy without ’’fine-tuning”. 
Phenomenological studies so based are referred to as ’’intermediate inflation” Il28l . 

During the big bang, the scalar held would emerge far from equilibrium, and, like 
the Bose gas illustrated in Fig.9, equilibrate by going through a period of quantum 
turbulence. The vortex reconnections that maintain the vortex tangle produce jets of 
energy that can create matter, like the solar flares created by magnetic reconnections 
illustrated in Fig. 3. This provides a new scenario for the inflation era, in which all 
matter in the present universe were created in quantum turbulence. 

A detailed analysis of the turbulent era based on the coupled Einstein-scalar equa¬ 
tions, unfortunately, faces the formidable problem of a non-uniform metric. We can 
bypass this difficulty with a phenomenological approach using Vinen’s equation, which 
deals with a uniform distribution of vortex lines. It is then possible to construct a sce¬ 
nario for the growth and decay of a vortex tangle, as illustrated in Fig. 14. Parameters 
can be so chosen that the lifetime of the vortex tangle is of the order of 10“^®s, during 
which time the universe expanded by a factor of 10^^, and a total amount of matter was 
created equal to what is presently observed, roughly 10^^ suns. This presents a new 
picture of inflation El- 

The big-bang era, including inflation and the creation of all matter through quantum 
turbulence, lasted about 10“^®s. After this, nonuniformities in the universe become 
important, and the Robertson-Walker metric ceases to be valid. The Halpern-Huang 


13 












potential, which is a high-cutoff approximation, may need corrections. Also, while 
there was only one scale in the big-bang ear, the Planck scale of 10^^ GeV, a new scale 
appears with the emergence of matter, the QCD scale of 1 GeV. Vortices created during 
the big bang era, including those in the vortex tangle, have core sizes of the order of the 
Planck length of lO^^^cm, and they will co-expand with the universe. Those created 
after the emergence of the matter, however, will have much smaller fixed core sizes that 
do not expand with the universe. 

After the big bang era, our model will be replaced by a cosmic perturbation the¬ 
ory, which must include the scalar field, and this will govern the creation of the CMB. 
The widely used KCDM model would be modified with inclusion of a vacuum com¬ 
plex scalar field, which supplies a cosmological-constant (A) via its energy-momentum 
tenor, and generates cold dark matter (CDM) through density variations of the scalar 
field. A complete reformulation is yet to be done, but we expect that the results of 
the usual cosmic perturbation theories will be largely preserved. The only new ele¬ 
ment would be quantum vorticity, which will make contributions to the tensor mode, 
in competition with gravitational waves. 

In the present universe, galaxies will attract superfluid to form halos of denser- 
than-vacuum superfluid around it, and these will be perceived as dark matter, through 
gravitational lensing n. When galaxy clusters collide, their haloes undergo distor¬ 
tions according to superfluid hydrodynamics, as observed in the ’’bullet cluster” Il29l . 
Computer simulations based on the NLKG may be found in lH. The literature contains 
studies of cosmic vorticity under the name ’’cosmic strings” ll^ |l33l. 

A fast-rotating body such as a black hole will drag the surrounding superfluid into 
rotation through the creation of vortices. The so-called ”non-thermal” filaments 041 
observed near the center of the Milky Way could be vortex lines, as shown in the upper 
panel of Fig.l5. The lower panels shows 3D simulations based on the NLKG ||3]. 

The core sizes of remnant vortices from quantum turbulence in the big bang era, 
originally of order of the Planck size of lO^^^cm, will co-expand with the universe, and 
in the 13.7 billion years since, they could reach sizes of the order of 10^ light years. 
Since matter was created in the scalar field, these cores are devoid of matter, and show 
up as voids in the galactic distribution. The so-called ’’stick man” configuration llT5l is 
shown in Fig.16, with a simulation by superposition of vortex tubes. 


7 Galaxy formation 

The CMB was formed about 10^ yrs after the big bang. Between that and the present, 
there was a long period of galactic formation, in which quantum vorticity may play an 
important role. A speculation of Lathrop ll3^ is that a galaxy can form from a large 
vortex ring with accretion of dust, which gravitates to form a central mass, squeezing 
the vortex ring into a shape with spiral arms, as illustrated in Fig. 17. Can one find any 
hint of this mechanism in currently-observed galactic properties? 

The central mass would correspond to the black hole observed at the center of all 
galaxies, whose mass M bears power-law relations to other galactic properties OTl : 

M (24) 
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Figure 15; Upper panel; The ”non-themal filaments” observed near the center of the 
Milky Way, Ref. IlT^ . could be vortex lines surrounding black holes. Lower panel; 3D 
simulations of a vortex-ring assembly arising from a rotating body at the center, based 
on the NLKG, in various perspectives. Ref. ||3l . 



Figure 16; Left; Voids in galaxy distribution in the ’’stick man” configuration. Ref. 
II 35 I . Right; Simulation by the superposition of three primordial vortex tubes, which 
have co-expanded with the universe to gigantic sizes. The inside of the tubes are devoid 
of galaxies, which cling to the outside dues to hydrodynamic pressure, like the metallic 
dust clinging to vortex tubes in superfluid helium shown in Fig.8. 
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Figure 17; Lathrop’s suggestion for galaxy formation; dust particles accrete onto a 
vortex ring, gravitate and clump, squeezing the ring into a spiral shape with a central 
mass, which would become a black hole. 


as indicated in the following; 


Stellar mass 
Luminosity 
Stellar velocity 


X p 
m 1.05 + 0.11 
L 1.11+0.13 
V 5.57 + 0.33 


(25) 


Assume that the dust particles initially accrete onto the vortex ring uniformly, and 
then clump up under gravitation. Assume further that a hxed fraction forms the central 
mass, which becomes a black hole and goes dark, while the rest remains luminous. 
This would mean m oc M, and L cc m, which are consistent with the relevant exponents 
being unity. 

The initial vortex ring may be generated by self-avoiding random walk (SAW) . 
(The fact that it is a close ring matters little for the arguments here.) Thus, the dust 
particles, which initially adhere uniformly to the vortex ring, form a SAW sequence, 
with the relation N ~ where R is the spatial extension, and N the number of steps. 
In our case, R corresponds to the size of the galaxy, and N is proportional to M + m, 
hence to M . This means that M scales like R^^^. Now assume that the total angular 
momentum J, which is a constant of the motion, scales as the galactic volume; 


M ~ R^'^ 

J ~ R^ (26) 

Dehning the stellar velocity v through J - Rmv oc RMv, we obtain 

M ~v^ (27) 

which is not inconsistent with observations. As an interesting note, the SAW exponent 
5/3 is the same as the Kolmogorov exponent in classical turbulence, and the Flory 
exponent for polymers ESI. 
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